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Abstract: We study perturbations around some cosmo logical backgrounds in the 
dRGT theory of massive gravity. We develop a general formalism to calculate the 
perturbations around any background. We derive the Lagrangian for fluctuations in 
the small scale limit, and for the open FRW solution we repeat the analysis around 
the full background. We find that the perturbations display similar properties: the 
longitudinal modes of the massive graviton are instantaneous at quadratic level, but 
they acquire a kinetic term at cubic order. 



Contents 

1 Introduction 2 

2 Action of Massive Gravity 4 

3 A formalism for perturbations 6 

4 Small scale expansion 7 

5 Cosmological solutions and perturbations 8 

5.1 de Sitter solution 9 

5.2 Flat FRW solution 10 

5.3 Open FRW solution 11 

5.4 Cosmological perturbations 13 



6 Conclusions 



14 



1 Introduction 

Recently, the possibility of constructing a consistent, Lorentz-invariant theory of mas- 
sive gravity has received renewed interest. Historically, the first quadratic Lagrangian 
which describes a massive spin-2 field was written by Fierz and Pauli [1]. Such a La- 
grangian propagates 5 degrees of freedom on a Minkowski background, although the 
mass term (in unitary gauge) breaks all four linearized diffeomorphisms which are a 
symmetry of the massless theory. However, Boulware and Deser (BD) [2] showed that 
the sixth mode, which is a ghost, reappears beyond the quadratic order. A renewed 
theoretical interest was spun by the Stueckelberg field formalism introduced in [3, 4], 
which allows a better understanding of the properties of the theory in terms of the 
effective theory of the helicity-1 and helicity-0 degrees of freedom. It becomes apparent 
that the BD ghost arises from the fact that the helicity-0 field satisfies a fourth-order 
equation of motion, and its absence in the quadratic Fierz-Pauli Lagrangian is guaran- 
teed by the "sick" kinetic term being a total derivative. One might therefore wonder 
whether a similar tuning can be found for the non-linear terms in the graviton La- 
grangian: the answer is yes, as was recently shown by de Rham and Gabadadze [5]. 
Finally, it was shown by de Rham, Gabadadze and Tolley (dRGT) that it is indeed 
possible to resum the whole perturbative Lagrangian into a 2-parameter theory, which 
is fully non-linear and Lorentz-invariant [6]. By construction, such a theory propagates 
only five modes around Minkowski background, and Hassan and Rosen [7] proved that 
the theory is indeed always free of the sixth mode, as the lapse function remains a La- 
grange multiplier after a redefinition of the shift vector. It is worth stressing that this 
does not guarantee that the theory is completely free of ghosts. In fact, around some 
backgrounds some of the five propagating degrees of freedom might become ghosts or 
show other instabilities. 

In dRGT theory, exact solutions with high degree of symmetry have been found: 
we refer the reader to [8-16] for some spherically symmetric and cosmological solutions. 
Some interesting features of these solutions are worth to be pointed out. For instance, 
spherically symmetric, static solutions in asymptotically flat space have a physical 
singularity at the horizon [10]. This is not true if one does not require asymptotic 
flatness [9, 14]. Perhaps, the most striking peculiarity of the dRGT massive gravity 
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is that there are no spatially flat homogeneous and isotropic cosmological solutions 
[11]: the same constraint that forbids the BD ghost does not allow for a non-trivial 
cosmological evolution of a flat FRW ansatz. Allowing for spatial curvature, it was 
realized in [12] that two branches of solutions to the constraint equation open up, both 
describing an open FRW metric. 

An interesting problem to investigate is the behavior of fluctuations around non- 
trivial backgrounds, because their stability is not guaranteed a priori. Since the po- 
tential of dRGT is written in terms of matrix square roots, a general perturbative 
expansion presents some technical difficulty. Up to now, perturbations have been stud- 
ied in the decoupling limit around Lorentz-invariant backgrounds [17]; the issue of 
superluminality has been considered in [18, 19]; linear cosmological perturbations have 
been studied in [20, 21] and an analysis of perturbations has been performed around 
the black hole solution in [14]. 

In this work, we take a small step to understand the behavior of the perturbations 
around some Lorentz-breaking backgrounds. First, we develop a general and algorith- 
mic formalism for writing down the fluctuation Lagrangian around generic backgrounds. 
We apply the formalism to the three cosmological solutions presented in [9, 11, 12]. We 
choose a local inertial coordinate system and we focus on high frequency, high momen- 
tum fluctuations in a small region of space-time. Our results are somewhat surprising. 
Around the three backgrounds, the quadratic Lagrangians for perturbations take essen- 
tially the same form, and the only propagating modes are the GR ones. The additional 
degrees of freedom, in fact, lose their time- kinetic term and their speed of sound is 
formally infinite. This degeneracy is broken at cubic order, because interactions induce 
a time-kinetic term for the additional modes. For the open FRW solution, we also 
consider the perturbative Lagrangian around the full non-linear background. In uni- 
tary gauge, we find that the quadratic Lagrangian only propagates the two helicity-2 
polarizations as in GR. At cubic level the degeneracy is broken and the 3 additional 
modes will propagate. 

A similar situation has been encountered in the decoupling limit lagrangian of 
massive gravity in [14, 17], in cosmological solutions of the bimetric massive gravity 
[22], and in the non-linear Fierz-Pauli theory in [23]. Also other modified gravity 
theories, for instance Hofava gravity [24], present similar problems. Recently, [21] 
showed how the vanishing quadratic time-kinetic term can be obtained as the isotropic 
limit of perturbations of an anisotropic Bianchi universe. They also showed that the 
modes propagate at the quadratic level when a no n- vanishing anisotropy is present, 
and one of them is always a ghost. 
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At classical level, the kinetic term can have either sign because it is a cubic term. 
But even if the square of the frequency has a negative sign, for wavelengths that are 
not very long the frequency of perturbations would be very large, and therefore above 
the cutoff of the effective theory. At the quantum level, a quadratic kinetic term can 
be generated, as there is no symmetry forbidding it, and its sign will depend on the 
interactions and parameters of the theory. However, we expect its coefficient to depend 
on the cutoff used to regulate the theory. Another issue to consider is the coupling to 
matter of such modes, and how do they modify the gravitational interactions among 
bodies. These are interesting and deep questions, and we will not address them in this 
paper, leaving their study for future work. 

The paper is organized as follows: in Section 2 we present the theory. In Section 
3 we develop the formalism we use to get the fluctuation Lagrangian, and in Section 4 
we explain the physical setup of our analysis. In Section 5 we describe the cosmological 
backgrounds we are interested in and derive the Lagrangians for perturbations. Finally, 
in Section 6, we discuss our results. 

2 Action of Massive Gravity 

In general, a Lorentz-invariant action for a theory of massive gravity can be written as 

/if 2 r 

S = ^J d4 *v^ [R - U(g^ r)} , (2-1) 

where U is a potential which depends on the metric g^ and on four auxiliary scalar 
fields <j) a , which are introduced in order to keep manifest diffeomorphism invariance 
[25]. The a 's are scalars under diffeomorphisms, but they transform as a Lorentz 
vector under an internal global SO(3, 1) symmetry group, <fi a (x) — > A a b <f) b (x). Matter 
fields are coupled to g^ u in unitary gauge, defined as (f> a = S^x* 1 : since the matter action 
is diff- invariant, the scalar fields are not directly coupled to matter. The most general 
action of the form (2.1) propagates two graviton polarizations plus four other degrees of 
freedom. Five of these correspond to the polarizations of a massive graviton, whereas 
the sixth mode is a ghost, known as the Boulware-Deser ghost [2]. However, recently 
there has been a proposal for a theory of massive gravity [5, 6], which has been shown 
to propagate just the five polarizations of a massive spin-2 field in [7, 26, 27] (for a 
review of massive gravity, see [28] ) . 

The non-linear potential presented in [6] describes a two-parameter family of the- 
ories, and can be written in the form 

U(g,„, r) = -m 2 [2e 2 (/C) + 6/3 3 e 3 (/C) + 24/3 4 e 4 (/C)] , (2.2) 
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where fC is the matrix defined by 

K% = 61 - VfT v , (2.3) 
and E is the following combination of the metric and scalar fields: 

W u = g» a d a <f> a d„<j> c r)ac- (2.4) 

Here i]^ denotes the Minkowski metric, which is taken as a fixed metric in the target 
space 1 . We denote by ej(/C) the following invariants of a 4 x 4 metric: 

ei(/C) = Tr/C (2.5a) 
e 2 (/C) = ^[(Tr/C) 2 -Tr/C 2 ] (2.5b) 

e 3 (/C) = - [(Tr/C) 3 - 3Tr/CTr/C 2 + 2Tr/C 3 ] (2.5c) 
e 4 (/C) = det/C, (2.5d) 

while e k (JC) = for k > 4. 

The expression (2.2) written in terms of the matrix /C is useful as /C starts linear 
in perturbations around the vacuum solution g^ u = r)^, <p a = For our purposes, 

we will find convenient to rewrite the action directly in terms of the matrix E^,: 

M 2 f - r r 1 i 

3 = -^ j d 4 x^ |i2 + 2m 2 \a + aiei(VE) + a 2 e 2 (VE) + a 3 e 3 (VE)J j , (2.6) 

where the ctj are not independent, being related to the fa as 

a = 6 + 12/? 3 + 12/3 4 , «i = -(3 + 9/3 3 + 12/3 4 ), 
a 2 = l + 6/3 3 + 12/3 4 , a 3 = -3(/3 3 + 4/3 4 ). 

Notice that in eq. (2.6) the e 4 term appears because we explicitly expand the invariants 
of fC. Also, we can safely drop the term e 4 (v / S) = Vdet E, since it is a total derivative 
and does not contribute to the equations of motion. 



One can modify the theory by considering a generic reference metric, or by considering the second 
metric to be dynamical as well, in which case the spectrum consists of a massless and a massive 
graviton [29]. 
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3 A formalism for perturbations 



The presence of the square root structure in the action (2.6) makes the analysis of the 
theory technically difficult. The problem of writing down the Lagrangian for perturba- 
tions around any given background is just a technical one: the expansion of the matrix 
square roots is difficult because the background matrix E and its perturbations 5E will 
not commute in general. One can rely on a Lagrange multiplier formulation of the 
action [6], or on the formulation in terms of vierbeine [30], but the calculations quickly 
become cumbersome because of the increased number of fields. Here we work out alge- 
braic equations which relate the invariants of the matrix \/E in terms of the invariants 
of E, which are readily perturbed (our method generalizes the calculations presented 
in [18]). These equations can then be solved order by order to get the Lagrangian for 
fluctuations. 

We start by establishing some notation. Denoting by Aj the eigenvalues of the 
matrix E, by definition the invariants e;(E) are given by 




(3.1a) 




(3.1b) 



i<j 




(3.1c) 



i<j<k 




(3.1d) 



Therefore, the invariants of the square root e^v^E) have the following expressions: 




(3.2a) 




(3.2b) 



i<j 




(3.2c) 



i<j<k 




(3.2d) 



It is straightforward to derive the relations: 




tt = 8i + 2t 2 , 



(3.3a) 
(3.3b) 
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Finally, one can combine the above identities to get the following quartic equations for 
the variables t 2 , t 2 , t\: 

t\ - A Sl t\ + (Qs\ - 8s 2 - 48v^I)^ 

+ 4(-s? + 4sis 2 - 16s 3 + 8s ly /sl)tl + (sj - 4s 2 + 8y/sl) 2 = , 
t\ - 2(s 2 + Qyfsl)tl - 8(s 3 + s 1 ^sl)t 2 + (s 2 - 2^) 2 - 4 Sl s 3 = , 
ijj - 4s 3 tl + {6s 2 3 - 8s 2 s 4 - 484 /2 )^ 

+ 4(-s^ + 4s 2 s 3 s 4 - 16s lS 2 4 + 8s 3 4 /2 )t 3 + (s 3 - 4s 2 s 4 + 8s 3 / 2 ) 2 

In principle, these equations can be even solved exactly to express the lagrangian as a 
scalar function of traces of E, but the resulting expression would be quite cumbersome 
and of very little use. Instead, we will solve eqs. (3.4a, 3.4b, 3.4c) perturbatively. Given 
an exact background solution g^, (fi a , we define the perturbations to the metric and 
scalar fields as g^ v = g^ u + h^ u and <p a = <p a + 7r a2 . We can construct the background 
matrix = g^ a d a cj) a dy(j) a and the perturbations = S — E. The invariants will be 
expanded as Si = Si + Ssi and U = ii + Sti, where Si = ej(E) and U = e^v^E). It is 
straightforward to solve eqs. (3.4a, 3.4b, 3.4c) for 5ti order by order as a function of 
h^, 7r a . Such a calculation method is algorithmic and it can be easily performed with 
a symbolic mathematics program (we used Mathematica in the following calculations). 

4 Small scale expansion 

The perturbative technique outlined in the previous section is very general and can be 
used to derive the fluctuation Lagrangian in different physical setups. It can also be 
applied to the case in which the reference metric is different from Minkowski or in the 
bimetric massive gravity theory [29], with minimal modifications. 

Here, we want to analyze the degrees of freedom and the stability of the theory 
around some cosmological backgrounds, by focusing on the Lagrangian of the graviton 
and Stueckelberg fields in a small region of space-time [25, 31]. 

The basic idea is that, for distances much smaller than the curvature radius, we can 
go to a locally inertial coordinate system such that the (background) metric, around a 
point x M = 0, is a small deformation of flat space-time: 

9tiu(x) = r]^ - -(R^vp + R^ ua )x a x 13 + 0(x z ) , (4.1) 

2 Anothcr possible definitions for the perturbations of the Stueckelberg is via the formula 4> a (t, x l ) = 
</> a (i+7T°, x 1 +tt 1 ), in which case the 7r a 's have simpler transformation properties under diffeomorphisms. 



(3.4a) 
(3.4b) 
(3.4c) 
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where R a Ruv ls the Riemann tensor of the background solution. In massive gravity, 
any coordinate transformation will modify the Stueckelbergs as well, so we will expand 
them in a Taylor series 3 : 

0a = A a x n + B a flv X ll X V + 0(X 3 ) . (4.2) 

The matrix A ^ is actually the vierbein of the unitary gauge metric at the point x = 
[26], and it will be proportional to <5° only around Minkowski space-time. In particular, 
we are spontaneously breaking the local Lorentz invariance even at the point x = 0, 
which will have important consequences at the level of perturbations. In the next 
section, we will work out the lagrangian for perturbations around this background, 
keeping only the terms which survive the x M — > limit in the equations of motion. 



5 Cosmological solutions and perturbations 

We are interested to study the behavior of perturbations around three background 
cosmological solutions: the de Sitter solution of ref. [9], the flat FRW solution described 
in the appendix of [11], and the open FRW solution found in [12]. It has to be noted that 
the solutions of [9] and [11] are isotropic but inhomogeneous, since it is impossible to 
find a flat FRW solution in which the scalar field sector is homogeneous and isotropic, 
as shown in [11]. However, if the graviton mass is small (m <C H), these solutions 
approximate very well the GR ones, because of a cosmological Vainshtein mechanism. 
If one allows for spatial curvature, an open FRW solution does exist [12], because this 
is invariant under a 5*0(3, 1) group, which is the diagonal subgroup of the space-time 
and of the internal 5*0(3, l)'s. 

The cosmological solutions admit a coordinate system in which the metric is of the 
FRW form. In isotropic coordinates, it can be written as: 

ds * = _ dT 2 + 2 (T) (i _g^ 4)2 ! (J.!) 

where K < is the curvature constant (we consider only spatially flat or open solu- 
tions). We go to a local inertial frame by performing the following non-linear transfor- 
mation to new coordinates t, x [32]: 



r = t--H(t)x 2 , y = —7-7- 
2 w y a(t) 



1 H -H-x 2 



(5.2) 



3 The Stueckelberg fields are derivatively coupled, so in principle we should keep the third order in 
the expansion. However, we arc interested to derive the equations of motion for perturbations in the 
x — > limit, in which case the third order does not contribute [26] 
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which brings the metric into the following form (up to second order in x): 

1 



ds' = - 



(H(t)+H 2 (t)) 



X' 



dt 2 + 



i--H 2 {t){i-n K {t))x' 



dx z 



(5.3) 



where H = a/a is the Hubble constant and Qk = —K/{a 2 H 2 ) is the usual curvature 
parameter. It will be useful to work with coordinates which are a small conformal defor- 
mation of Minkowski space-time. In such a form the metric will deviate from Minkowski 
both for distances and times ~ H' 1 , and the gauge freedom will be completely fixed. 
The needed infinitesimal diffeomorphisms are [32]: 



£o = ~(2H + H 2 + n K H 2 ) Ux 2 + 



f tx * + ^ , £ = i(2if + H 2 + n K H 2 )xit 2 , (5.4) 



and the metric reads 
ds 2 



i - \{i - n K )H 2 x 2 + i ( l + n K + 2^ ) hH 2 



-dt 2 + df 5 



(5.5) 



Notice that in such a form, the parameters H, H and Qk are all evaluated at a fixed 
time. The different solutions will differ by the spatial curvature and by the Stueckelberg 
sectors, which we will now describe in detail. 



5.1 de Sitter solution 

We consider the solution discussed in the main section of [9], which is restricted to 
the one-parameter choice fa = — 4(3^. For notational simplicity, we will denote C = 
(1 — 12/?4)/(l — 8/S4) in the following formulae. In FRW coordinates, the metric is de 
Sitter in planar coordinates, while the Stueckelbergs are given by 

~l + H 2 p 2 -e- 2HT ~ 



2C 
3H5 



arctanh 



1 - H 2 p 2 + e~ 2HT 



-p[ W e \y\ 



^f"V, ( 5 - 6 ) 



where 5 is a dimensionless integration constant in the range < 5 < 4C 2 /9, and the 
function p(z) is defined by the equation 



2 ^ TT V16C 4 - 81<5 2 + 36C 2 5 2 H 2 z 2 
P {Z) = -36 CHZ W - 9 



(5.7) 



Using the transformation to inertial coordinates (5.2) and (5.4), the Stueckelberg back- 
ground up to second order in x becomes, up to an irrelevant constant, 

1/2 „ 

x 2 , = —x' . (5.8) 



2C 3H /16C 4 

-t+—[ " 1 



35 AC V 815 2 



2C 
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The quadratic Lagrangian for the perturbations reads: 



m 2 M 2 l 
L-1 — 



+ 



m 2 Ml x 



12^ 
- - 12/3 4 



4C 2 
~9T 



-F 2 



2C 



hijidiTTj - dkiikSij) + 



4 



/i 2 



+ 



4C 2 - 96 
AC 2 + 95 



AC 



Hx l F-7r- 



(5.9) 

where we have defined = d(Kj — <9j7Tj. 

In the high momentum limit, the mixing with gravity can be ignored. In the 
Stueckelberg Lagrangian, it is clear that none of the 7r's is propagating, their equations 
of motion involving only spatial derivatives. A particularly bad situation arises in the 
Lorentz-invariant limit, i.e. when the parameter 5 — > 4C 2 /9, since in this limit the full 
second-order Lagrangian vanishes: here perturbations are infinitely strongly coupled. 

The cubic lagrangian reads: 

(1 - 12/3 4 )C 



545 



96 F^ 



(AC 2 + 95)F ij F jk d (i TT k) - 365FijdiK\ 



dk^k) + 3245 2 (<9 (i 7r j) - d h -K k 5 ij )(i: i 'kj 



27T i <9 ? -7r° 



(5.10) 



from which it is clear that there is a time- kinetic term for the 7Tj, while 7r° is algebraically 
determined. One might think that this result is an artifact of the choice of gauge, but it 
is not. In fact, at quadratic level we wrote down the full Lagrangian for perturbations. 
Even if we choose a gauge in which the perturbations are only in h^, it is clear that 
h 00 and h 0i do not appear and the constraint structure is the same as GR. At cubic 
order, we checked that h 00 appears as a Lagrange multiplier in the mass term, while 
hoi appear quadratically, thus modifying the GR constraints and giving a kinetic term 
to three additional modes. Another doubt one may have is that the cubic order kinetic 
term could be removed by a field redefinition (we refer the reader to [33] for a discussion 
on this point). By writing down the equations of motion, it appears that two initial 
data are required to specify the solution for 7Tj. Therefore, we conclude that three 
longitudinal modes are propagating at cubic order. 

5.2 Flat FRW solution 

A generic flat FRW solution has been derived in the Appendix of [11], for the parameter 
choice /?3 = /?4 = 0. The scale factor is determined by the matter content of the 
Universe, and the Stueckelberg sector is, in FRW coordinates, 

9 r ds , s „ A 9f \ ,, 3 



16T 



a(s)H(s 



+ cl(t)T 1 



16T 2 



4> l = W 



(5.11) 
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where T is an integration constant with dimensions of time 4 . Using the transformation 
to local inertial frame, up to 0(H 2 x 2 ), the scalar background becomes 



9 + 16a 2 H 2 T 2 



-t 



16a 2 H 2 T 2 - 9 



t 2 T 

{H + H/H)--H X - 



16aHT 16aHT 
The perturbative Lagrangian at quadratic level reads 



■x\ (5.12) 



m 2 M Pl (3 - AaHTf 
' 2 24aHT 
m 2 Ml x 3 - AaHT 
+ 2 3 + AaHT 



13 9 

1 F ij ~ 2 h ij( d iKj ~ SijdkTtk) + ( h ij ~ h tk) 



+ 



1 6n 2 H 2 T 2 4- 9 

AH^d 1 Tl j - Hx'FijdjTT ; TT J Hx'FaTTi 

03 24aHT 3 3 



(5.13) 



which has the very same structure of (5.9), and does not propagate any degrees of 
freedom. 

At cubic order, the Lagrangian reads 
A = (3 ^f Q 5P FjI'jiA,-!.-, + h^ih^h, - \f 2 (vr° + d k TTk) 



lUaHT 
+ 



lQaHT 



3 J J 6 
(d(i7Tj) - d k 7i k Sij) (itiitj + 27r;<9j7r° + di^djir ) , (5.14) 



(3 + AaHT) 2 

which is very similar to (5.10), and the same considerations apply. 
5.3 Open FRW solution 

For the solution presented in [12], the scalar field background reads (in FRW isotropic 
coordinates) 



c±a(T) 



y 2 



c±a{T) 



\K\ (1 - \K\f/A) 2 ' r 
There are two branches of solutions identified by the constant c±: 

7i + 72 ± v 7 7? - 72 



K|y2/4 



(5.15) 



72 



where, to simplify the following expressions, it is useful to define: 

7l = l + 3/3 3 , 72 = 3(/3 3 + 4/3 4 ). 



(5.16) 



(5.17) 



4 This is one particular solution of a non-linear differential equation, every solution of which gives 
an equally valid Stucckelberg sector for the FRW metric. For instance, taking 04 = 0, the dS solution 
of the previous section obeys such equation. 
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Since the theory is defined by the positive square roots, we require c± > 0. In the local 
inertial frame, up to second order in Hx, we get 



c± 1 c+ 
7==* + o^= 



(H + ij/i/)t 2 - (1 - n K )Hx 2 



c±x 



(5.18) 



As before, we write the quadratic Lagrangian, keeping only the terms which survive 
the x — > limit in the equations of motion: 



c _ m 2 Ml x (1 



K 



4c± 



F„ - h i:j (diiTj - d k n k 5ij) + —(h^ - h 



+ 



+ 



m 2 M^ (1 - VQ^)A 



-AH^djTTj + Hx^ijdjTi + -^=Hx i F ij 7i j 



A 



(5.19) 



where we have defined A = (1 + 271 + 72 — c±(7i + 72)). Again, this Lagrangian has a 
similar structure to (5.9) and (5.13) and describes just non-propagating modes. This 
result resonates with the conclusions of [20], who analyzed the cosmological pertur- 
bations at linear order and found instantaneous propagation of the vector and scalar 
modes. We will extend their analysis around the fully non-linear background in the 
next section. 

An interesting observation can be made regarding the Qk — > 1 limit. In fact, 
the choice Qk — 1 corresponds to the Milne universe solution, which has the metric 
ds 2 = — dr 2 +r 2 dx 2 in FRW coordinates, but it is just a reparametrization of Minkowski 
space-time (in fact, the metric becomes exact Minkowski in local coordinates). The 
quadratic Lagrangian (5.19) in this limit turns out to be C = 0, which signals that the 
perturbations become infinitely strongly coupled. This has to be compared with the 
healthy lagrangian £ = —\F 2 U + . . . around the Minkowski solution. The difference in 
behavior should not surprise us. In fact, the scalar equations of motion for the open 
FRW ansatz have three branches of solutions [12]: two are present only in the case 
K 7^ and correspond to the non-trivial cosmological evolution we just described, and 
the third one forbids any evolution of the scale factor, so that the metric reduces to 
flat space-time. So, perturbing around the "cosmological" branch and then taking the 
Minkowski limit does not commute with considering the Minkowski branch and then 
perturbing around it: the first choice displays a bad behavior. 
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At cubic order, the Lagrangian reads 
£ 3 = ~ — "7^— [(7i + 72 - 3c ± 7 2 )F ij F ifc 9( fc vr i) - c^FijdkTTk] 

- g [48 (d m - d k -K k 5 l3 ) KidjTr - 3 (G% - G 2 kk ) tt° + 2G ij G jk G ki - ?,G%G kk + G 2 kk ] 
7i + 72 - 2c±7 2 ( „ a o . 1 2 . o 



+ 71 + 72 — 2C±72 a/Ur- - d k Ti k 5ij) (wi-kj + 2^^-77° + 0i7r°9j7r°) , 

(5.20) 

where Gjj = 2d^7Tj), and for simplicity we performed an expansion in the small pa- 
rameter VLk-i retaining the leading order for each term. Again, the it is algebraically 
determined while the 7Tj acquire a time-kinetic term. 

5.4 Cosmological perturbations 

The analysis of perturbations on small scales highlights the degeneracy of the cosmo- 
logical backgrounds. To convince ourselves that this is not an artifact of the approx- 
imations made, it is also worth considering the behavior of perturbations also around 
the fully non-linear background solutions, without any small-scale expansion. This is 
easy to do in the case of the open FRW solution, which we will now study in detail. 
We start by writing the background in the following coordinate system: 

ds 2 = -dt 2 + a 2 (t)n ij dx l dx j = -dt 2 + a *® (tfSij - ./•'./-' )d/-'(l.r' , (5.21) 

where, K < is the value of the spatial curvature. The background Stueckelbergs are 
of the form 

10 c±a(t) 



a/1 + \K\x*, ft = c±a(*y . (5.22) 



We then choose the unitary gauge for the perturbations, i.e. we consider g^ u = gj/ + 
h^y and we leave the a 's unperturbed. Applying the method of section 3, we find the 
mass term of the Lagrangian at quadratic order to be 



13 n 



K 



c± [(71 + 7 2 )c± - 3(1 + 2 7 i + 72)] {h 2 - h tj h tj ) (5.23) 
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where we have defined h = Q^hij, = Cl tk QP hki- The term C c , c . is the effective 
cosmological term —M^Aesy/^g expanded at second order, where 

o 

Acff = [2c ± ( 7l 2 - 72) - (7i + H ~ 272)] • (5.24) 

72 

From this expression it is clear that in the full Lagrangian only two degrees of freedom 
are propagating, as in GR. In fact, Hqq and h 0i do not show up at all in the mass term, 
so that the full Lagrangian has the same constraint structure as the GR one. We can 
now ask if this enhanced symmetry is respected by the cubic order Lagrangian, and 
the answer is no. In fact, the expression for the mass term at cubic order reads 



M^m 2 a 3 Vdettt 
1672v / n^ 



± (67? + 7 2 2 - 572 + 27172) - (7i + 7 2 )(1 + 2 7 i + 72)) hi^hi 



3 



+ ((7i + 72)(1 + 27i + 72) - 2c±( 7l 2 - 72 )) h^h - ^ 72 (1 + 3 7l + l2 )h 3 
+ ((71 + 72XI + 271 + 72) - c±(2 7l 2 + 2 7l72 - 72 + 72 2 )) hoo(h ij h i:i - h 2 " 

+ An K (( 7l + 72 )(1 + 2 7l + 72 ) - c ± (2 7l 2 + 2 7l72 - 72 + 72 2 )) (h ij - hQ ij )h 0i h 0j 

(5.25) 

where, for simplicity, we only display the leading terms in an expansion in Vtx for each 
coefficient. Looking at this expression, we see that, while h 00 appears as a Lagrange 
multiplier, h 0i appears quadratically, which means that the cubic Lagrangian propa- 
gates three more modes than the quadratic one. This is in fact analogous to what we 
have seen in the small scale limit 5 . 



6 Conclusions 

The cosmological solutions we have considered look different from each other. However, 
they share similar properties at the level of perturbations. Namely, longitudinal modes 
of the graviton do not propagate at quadratic order in perturbations, but they do get 
a time-kinetic term by cubic order interactions. 

This means that, schematically, their dispersion relation is of the form k 2 = at 
quadratic level, but it is modified to k 2 = eu 2 when considering the cubic Lagrangian. 
Here e is actually the fluctuation field, whose sign in particular is ambiguous. There- 
fore, perturbations around the cosmological background can have negative frequency 

5 Notice that the term involving /loi^oj appears with a coefficient proportional to s/TTk, as the term 
involving ■ki'kj in eq. (5.20). This is a consistency check for our calculations. 
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squared, but for short wavelength the frequency is actually very high, and above the 
cutoff of the effective theory. Quantum effects can generate a quadratic kinetic term, 
whose sign will depend on the parameters and interactions of the theory. We expect 
that its coefficient will depend on the cutoff used to regulate the theory. An important 
question to address is to include the coupling to matter, and to understand the exper- 
imental consequences of the behavior of the longitudinal modes. In general, it would 
be interesting to understand what is the reason for the vanishing of the kinetic terms, 
and in particular how this behavior is related to the symmetries of the auxiliary and 
physical metric. 

All these questions are important and interesting, and we leave them for future 
studies. 
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